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ABSTRACT 

r^ I The S stars orbiting the Galactic center black hole reach speeds of up to a few percent the speed of 

^ ■ light during pericenter passage. This makes, for example, S2 at pericenter much more relativistic than 

Oj. known binary pulsars, and opens up new possibilities for testing general relativity. This paper develops 

' ' a technique for fitting nearly-Keplerian orbits with perturbations from Schwarzschild curvature, frame 

^ • dragging, and the black hole spin-induced quadrupole moment, to redshift measurements distributed 

c/3 , along the orbit but concentrated around pericenter. Both orbital and light-path effects are taken into 

account. It turns out that absolute calibration of rest-frame frequency is not required. Hence, if pulsars 

on orbits similar to the S stars are discovered, the technique described here can be applied without 

V^J I change, allowing the much greater accuracies of pulsar timing to be taken advantage of. For example, 

I — ' pulse timing of 3 fis over one hour amounts to an effective redshift precision of 30 cm s^^, enough to 

^^ , measure frame dragging and the quadrupole moment from an S2-like orbit, provided problems like the 

^^ ' Newtonian "foreground" due to other masses can be overcome. On the other hand, if stars with orbital 

O ' . . 

• ' periods of order a month are discovered, the same could be accomplished with stellar spectroscopy from 

^ ! the E-ELT at the level of 1 km s"^ 

O 



C^ 



X 



Introduction 



Tracing the orbits of th e S stars in the Galactic center region reveals a central mass of ~ 4 x IO^Mq (jGhez et al. 



Ci I l2008t iGillessen et al.ll2009a ). presumably a supermassive black hole. The roughly 20 known S stars are mostly main 



sequence B stars, and it is likely that they share their environment with many fainter stars. Hence discovery by 
the next generation of telescopes of stars far closer to the black hole is anticipated. The dynamics of the orbit and 
light trajectories provide an opportunity to test hitherto unobserved predictions of general relativi ty in the black 



hole neighborhood, and in doing so, confirm or revise our understanding of the nature of gravity ( Gillessen et al. 



20101) . Meanwhile, searches for pulsars in the Galactic-center region are also underway (e.g.. lMacauart et al.ll2010 ) 



If pulsars on orbits similar to S stars are discovered, their time-of-arrival measurements may provide accuracies 
orders of magnitude beyond those available from optical spectroscopy. 



It is interesting to compare S stars with other relativistic orbits. General relativistic perturbations of a nearly 
Keplerian orbit typically scale as an inverse power of the classical angular momentum ■\/a(l — e^). So a simple 
measure of the strength of relativistic effects in a nearly Keplerian orbit is the dimensionless pericenter distanceo 
Figure [1] plots pericenter distance against orbital period for a sample of S stars and binary pulsars, together with 
Mercury and artificial satellit es. At p resen t, binary pulsars are the leading laboratories for general relativity. The 
principles are summarized in iTavlorl (|1994I ). Taylor's Figure 11, depicting five intersecting curves, indicating four 
distinct tests of general relati vity, is especi ally memorable, and descendants of this figure with newer data appear in 
recent work (e.g., Figure 3 in lKramer fc Wex 20091). So it is remarkable that the star S2 has a pericenter distance 



~ 3 X 10'^, an order of magnitude lower than the most relativistic binary-pulsar systems, and four orders of magnitude 
lower than Mercury. This suggests that Galactic center orbits may show relativistic effects not measured in binary 
pulsars, such as frame dragging. On the other hand, 82 has an orbital period of ~ 15.9 yr versus one day or less 
for some binary pulsars. Thus, for the known S stars, or pulsars on similar orbits, we cannot build up a relativistic 
signal over many orbits. A different strategy is needed, based on detailed observation of one or a few orbits, and 
paying special attention to pericenter passage, when relativistic effects are strongest. 

Several different possibilities for observing relativistic effects have been discussed in the literature. To see how 
these relate to each other, let us consider the different physical effects involved in observing a relativistic orbit over 
time. These are illustrated in Figure [2] Shown schematically in this figure is a star on a precessing orbit, which 
emits photons (pulses, wavecrests) at a fixed frequency in proper time. These photons then take a curved path to 
the observer. When they arrive, both their frequency and their travel direction have been altered. What can the 
observer detecting these photons infer about the metric? This observer must take multiple physical processes into 
account. 



1. 



Time dilation at the source is a consequence of the equivalence principle, and is the strongest relativistic 
pertu rbation. Its measurability in the context of the Global Positioning System (GPS) is well known (jAshby 
20031 ). The manifestation of time dilation in pulsar timing is known as the Einstein time delay. For th e S 



stars, a spectroscopic detection of time dilation is expected around pericenter passage (jZucker et al.ll2006l ). 



2. Orbit perturbations are a test of space curvature, and at higher order, of black hole spin. 

• The leading order precession of orbital periapse is well measure d in binary pulsars. (For a discussion of 
how pulsar timing effects scale to the Galactic-center region, see lPfahl fc Loebll2004l . .) Astrometric mea- 



surement of precession in the orbit of S2 is considered feasible with current instruments (JEisenhauer et al 



20091 ) ■ but the long orbital period, and the likelihood of Newtonian precession due to the distributed mass. 



both present serious difficulties. If stars much further in are discovered , however, precession of orbita l 
planes due to higher-order spin and quadrupole effects become accessible ( Willll2008tlMerritt et al.ll2010l ). 



Orbital decay due to gravitational radiation is not considered measurable because of the long time scales 
and the extreme mass ratio between the star or pulsar and the supermassive black hole. 

In contrast to t he above secular effect s , there are also GR induced velocity perturbations, which vary 
along the orbit ( Kannan fc Sahall2009t IPreto fc Sahal 120091) . In binary pulsars, velocity perturbations 



*In this paper the semi-major axis a and the pericenter distance a(l — e) are always expressed in units of GM/c^ and are therefore 
dimensionless. 



are not considered especially interesting and are subsumed within Roemer delays. For Galactic center 
stars however, the situation is different. First, the velocity perturbations arc larger, and second, they 
are concentrated around pericenter passage, thus offering a better prospect for isolating relativistic from 
Newtonian perturbations. 

3. Light paths and travel times are affected by the black hole mass, and again at higher order, by spin. 

• Strong deflection of starlight ( Bozza &: Mancinil 12009 ) or a pulsar beam ( Wang et al.ll2009bl la[) would be 



a spectacular though rare event. 

Small perturbations of photon trajectories, too small to measure astrometrically, can nevertheless produce 
detectable changes in the light travel time. In the pulsar literature this is known as the Shapir o delay. 
The analogous effect for S stars is a time-dependent redshift contribution ([Angelil fc Sahall20I0l ). which 
for a star like S2 is comparable to the velocity perturbations. Hence, in the Galactic center, relativistic 
perturbations on both orbits and light paths must be considered. 



Later in this paper, we will group time dilation with the orbital effects. 

Motivated by the above considerations, in this paper we develop a method to fit a relativistic model to ob- 
servables. This strategy (i) takes relativistic perturbations on both the orbits and the light into account, (ii) treats 
stellar spectroscopy and pulse timing in a unified way, and (iii) is well suited to analyzing data obtained from a 
single orbital period or even less. We include the effects of time dilation, space curvature, frame dragging, and 
torquing-like effects induced by the quadrupolc moment from the black hole spin. Some further issues remain; most 
importantly, how to include Newtonian perturbations (from mass other than the black hole's) in the fit — hence 
the "towards" in the title. 

The observable we consider is the apparent frequency v. As remarked above, this can refer to spectral lines or 
pulses. The source frequency :^o is in principle known for stars, but unknown for pulsars. Note that 



cln(i'o/'^) = cln(l + z) ~ cz 



(1) 



where z is redshift in the usual definition. For this paper, however, we will use "redshift" to mean cln(^'o/i^). The 
possibly unknown source frequency now appears as a harmless additive co nstant in the redshift. For the S stars, 
the current spectroscopic accuracy is « 10 km s~^ under optimal conditions (JGillessen et al.ll2009a[ ). If a pulsar on a 
similar orbit is discovered, the accuracy would like ly be much higher. Taking, as an example, a one- hour observation 
with pulses timed to 3/iS (cf. IJanssen et al.l 120101 ) implies an accuracy of one part in 10^, equivalent to a redshift 
accuracy of 30 cm s^^. The same level of ac curacy is not inco nceivable from spectroscopy of S stars, since planet 
searches routinely achieve Af < 1 m s^^ (e.g. iLovis et al.ll2006| ). but would require some technical breakthroughs to 
achieve for faint infrared sources like the S stars. 

For a given orbit and redshift accuracy, the principal quantities we will calculate are the signal-to-noise ratio 
S/N for different relativistic effects. For each of four relativistic effects — time dilation, space curvature, frame 
dragging, and quadrupolc moment effects — we will present the results in two ways: (i) the redshift accuracy needed 
to reach some S/N for a given orbit, and (ii) the orbital parameters needed to reach some S/N for a given redshift 
accuracy. 



2. Model 

The main calculations in this paper will assume that the trajectories of both S stars (or pulsars) and photons 
are geodesic in a pure Kerr metric. Naturally, geodesies of the former arc timclike, and the latter null. Even if 
Einstein gravity is correct, the pure Kerr metric is naturally an approximation because this solution to the field 
equations neglects all mass in the vicinity of the black hole (more on this later) , and the mass of the star itself. In 
general relativity, geodesic motion is fully described by the super Hamiltonian 

H=\g^^''p^p,, (2) 

with H = {] for null geodesies. 

Rather than carry out all computations with the Hamiltonian resulting from the full Kerr metric, it is useful 
to consider two different approximations for the cases of orbits and light paths. These we will call iJstar and i7""". 
(Word labels in subscripts indicate orbits, in superscripts light paths.) We express these approximate Hamiltonians 
perturbatively as 

Hstar = -ffstatic + £ ^Kop + £ -^Schw + £ ^FD + £ ^q, (3) 

and 

where the meanings of the various component Hamiltonians will be explained shortly. Note that e is just a label to 
indicate orders: if e" appears, the term is of order w" (where v is the stellar velocity in light units), but numerically 
e = 1. The series in ([3|) and (|4]) actually consist of the same terms, but because if star is to be applied only to 
trajectories that are strongly timelike, and iJ""'^ only to compute null geodesies, the orders of same terms are often 
different. 

Expressions for all the Hamiltonian t erms are derived perturbatively from the Kerr metric, using the basic 
method explained in lAngelil fc Sahal ( 2010I ). We will not repeat the details here, but simply list the terms and the 



ir 



physical meanings. 

First we consider the orbit (Equation [3]). 

• To leading order, the star feels no gravity and the clock simply ticks: 



i^static — — ;r- (5) 



• 



At next-to-leading order, gravity first appears, with 

rr _Pr , pI , pI Pt (f.^ 

ffKcp-y + ^ + ^^^— ^--, (6) 

which is the classical Hamiltonian modified by 1/r — > p^/r. This modification leaves the problem spatially 
unchanged, but introduces a temporal stretching causing a gravitational time dilation. This is originally a 
consequence of the Einstein Equivalence Principle, and affects the redshift at 0(w^). 



• Appearing next are the leading order Scliwarzschild contributions, 

-nschw — 2 ' ^'' 

of which the first term effects a further time dilation, and the second curves space, leading to orbital precession. 

• Next, the first spin term emerges with 

^™ = -^. (8) 

producing frame dragging, an r-dependent precession around the spin axis. Here s is the spin parameter, 
which points in the 9 = Q oi +z direction. Maximal spin is s = 1. 

• The last set of orbital terms we will consider are 

which are the leading-order quadrupole moment terms. These result in the orbit being torqued towards the 
spin-plane, as well as other less intuitive effects. 

Then we consider light paths (Equation |4]). 

• At leading order, we have 

rrMink _ Pt , Pr , Po P4> (^r.\ 

which is to say, the spacetime is Minkowski and the photon trajectories are straight lines. 

• The leading order Schwarzschild contribution 

2 2 

rrSchw __^t Pr_ /ii\ 

r r 
(which notice is not identical to -ffsdiw) introduces lensing and time delays. 

• The frame-dragging term for photons 

debuts at one higher order than TJpD- Frame dragging lenses photons around the black hole in a twisted 
fashion. A 0-dependent time-dilation also occurs. 



• 



Finally, we consider higher order Schwarzschild terms, and quadrupole moment terms 

„q _ 2p^ s^ sin" Opl s^ cos^ dpi s^pl 

r2 2r2 2r4 2r4 sin" 6) ^ ' 

which, as before, produce a torque towards the spin plane and as well as some more subtle effects. 



The well-known separability properties of the Kerr metric (cf. IChandrasekhaiJ 119831 ) provide solutions up to 



quadratures, but not explicit solutions. Analytic solutions for geod esies and null geodesies are available for various 



special cases, for example for the leading-order Schwarzschild case (JD'Orazio fc Sahall201Cll ). But the present work 
uses numerical integration of Hamilton's equations. 



3. Calculating redshifts 



Given a set of orbital parameters for the star, wc calculate a redshift curve, meaning ln(i'o/;/) against observer 
time or pulse arrival time, by the following method. 

First the orbit is calculated by numerically solving Hamilton's equations for i/gtar (Equation [3]). If desired, 
particular orbital relativistic effects can be isolated by omitting other relativistic terms from i/gtar • The independent 
variable is, of course, not time but the affine parameter (say A). Along this orbit, from pairs of points separated by 
AA, photons are sent to the observer. The difference in proper time between the emission of two photons is 

Ar= ^ ^^ . (14) 

The photons themselves travel along paths determined by iJ""" (Equation |4|) with the additional condition iJ""!' = 
0. Again, one is free to isolate particular relativistic effects on the photons by discarding terms from i7""''. If Ai 
is the difference in the photons' times of arrival at the observer, then 

t/ At ^ ^ 

This calculation is carried out along the orbital path. Some snapshots of this procedure arc illustrated in Figure [2j 



Now, photons are emitted from the star in all directions, but we need to find exactly those photons which reach 
the observer. Solv ing this boundary valu e problem is the computationally intensive part. Our algorithm for doing 
so is explained in lAngelil fc Sahal (|2010l ). In stronger fields, the redshift curve takes longer to compute than in 
weaker ones, both because the orbit integration needs smaller stepsizes over more relativistic regions, and because 
the boundary value problem requires more iterations before satisfactory convergence is reached. 

The observable redshift is of course a combination of all the relativistic effects atop the classical redshift. One 
can, however, estimate the strengths of each effect in isolation, by toggling each term in Equations ^ and @ on 
and off, and then taking the difference in redshift. Figures [3] and |4] show relativistic redshift contributions measured 
in this way. These orbits refer to the pericenters of a range of orbits, varying in a but with e and the other orbital 
parameters fixed at the values of S2. These effects have simple scalings with the orbital period, readily deduced 
from t he scaling properties of the Hamiltonians, and summarized in Table [1] Numerical results in lAngelil fc Sahal 
( 20101 ) verify that such scalings apply not only at pericenter but all along the orbit. 



The extended mass distribution in the Galactic center region, due to all the other stars, stellar remnants and 
dark matter particles that are present, introduces Newtoni an perturbations. Th e distribution and normalization of 



this mass are poorly c onstrained on the scales of interest (jSchodel et al 



( Gillessen et al.ll2009bl ) and some numerical experiments have been done 
symmetry for this distributed mass, models of the form 

p (r) (X r^'' 



20091) ■ although so me upper limits exist 



Merritt et al.ll2010l ). Assuming spherical 



(16) 



are sometimes adopted. Such a model ignores the torques due a nonspherical or discrete mass distribution, and the 
only change it implies in the orbital dynamics is an additional (prograde) pericenter advance. Figure |5| compares 
the orbital and photon signals strengths for each effect, including an extended mass distribution of the type p6p . 
Note however, that in the analysis which follows, we do not attempt to include the effects of the extended mass 
distribution. 



4. Fitting 

The computational de mands for the presen t work are far greater than calculations of a few redshift curves 
for given parameters (as in lAngelil fc Sahall201(J ). because fitting redshift curves in a multidimensional parameter 
space requires evaluations of up to tens of thousands of such curves. To this end, parallel functionality was added 
to the implementation. The work is distributed evenly among the available processors. A workload is the charge 
of calculating the redshift from a single position on the star's orbit, i.e., a single evaluation of ([T5|). calling for the 
finding of two photons travelling from star to observer. The program is available as an online supplement. 



4.1. The parameters 

The redshift curve of an S star or pulsar depends on eight essential parameters, but there can be any number 
of additional parameters for secondary effects. In this paper we consider nine parameters in all, as follows. 

• The black hole mass Mbh sets the overall time scale, and accordingly we express it as a time. Changing A/bh 
simply stretches or shrinks the redshift curve in the horizontal direction. 

• The intrinsic frequency i^o for pulsars is the pulse frequency in proper time, whereas for spectroscopy, it has 
the interpretation of absolute calibration. Altering vq simply shifts the redshift curve vertically. 

• Then there are the Keplerian elements, referring to the instantaneous pure Keplerian orbit with the initial 
coordinates and momenta. Our orbit integrations all start at apocenter. 

1. The semi-major axis a (or equivalently the period P = 27ra"^'^A/BH) is single most important parameter 
dictating the strength of the relativistic signals. 

2. The eccentricity e sets how strongly peaked the redshift curve is at pericenter. 

3. The argument of pericenter w sets the level of asymmetry of the redshift curve. 

4. The orbital inclination / changes the amplitude of the redshift. Classically, the inclination enter the red- 
shift depends only on Mbh sin'^ / (on ^^InZ-^^totai ^'^^^^ ^ ^'^^' finite mass ratio). In relativity, the degeneracy 
is broken, because time dilation is independent of /. 

5. The epoch eo basically the zero of the observer's clock. 

The longitude of the ascending node il does not appear, because, in our chosen coordinate system, ft rotates 
the whole system about the line of sight, which has no effect on the redshift. Formally, we simply fix J7 = 0. 

• For the pulsar case, we include a simple spin-down model, with a constant spin-down rate vq ^^ v^ — v{ta). 
Here v is an additional parameter. 

We take the black hole spin as maximal, and to point perpendicular to the line of sight. Proper motion of the 
black hole is not considered. 

To fit, we use a quasi-Newton l imited- memory Broyden-Fletcher-Goldfarb-Shanno optimization routine with 



bounds (|Zhu et al.lll997l:lBvrd et al.lll994ll . 



4.2. Signal to noise 

With an orbit fitting algoritlini in hand, we now need to quantify the notion of detectability of particular 
relativistic effects. One way would be to attach a coefficient to each term in the relativistic Hamiltonians and then 
see how accurately that coefficient can be recovered. For this paper, however, we take a simpler approach. Since, 
for now, we only aim to identify which effects and which regimes are promising, we will simply attempt to estimate 
the threshold for detecting the presence of each term in the Hamiltonians ([3]) and @ . Accordingly, we proceed as 
follows. 

1. For some chosen parameters, we generate a redshift curve including all Hamiltonian terms up to some 0(e"). 
These redshift curves are sampled at 200 points, with dense sampling near pericenter. 

2. We add Gaussian noise to the redshifts, at some chosen level. 

3. We then fit the parameters, with a redshift curve lacking a particular Hamiltonian contribution. If the noise 
level is too high, a reduced x^ — 1 will be obtained, otherwise Xred ^^^^ ^e higher. We define 

S/N^^, (17) 

as the signal to noise. 

Naturally, it is necessary to check that a large Xred ^^ ^'^^ ^^e result of some algorithmic problem, by verifying that 
a good fit is obtained if the Hamiltonian term in question is consistently included. 



4.3. Redshift-accuracy demands for S2 

Figure IHl shows the S/N of time dilation, space curvature, frame dragging, and quadrupole terms, all as a 
function of redshift accuracy, assuming the orbit of S2. Time dilation is well above the current spectroscopic 
threshold of w 10 km s~^. Space curvature is somewhat below, while frame dragging and quadrupole are far below. 

The detectability thresholds are summarized in the first part of Table [21 



4.4. Orbit-size demands for given redshift accuracy 

In Figures [3 [8] and [9l the redshift accuracy is set to lOkms"^, Ikms^^ and 30cms~^ respectively, while the 
orbital period varies along the horizontal axis. 

Table [D summarizes the detectability thresholds. 



5. Summary and Outlook 

This paper addresses the problem shown schematically in Figure [21 which is to infer orbital parameters and 
detect relativistic effects from redshifts or pulsar timings along an orbit. The specific observable of interest is the 



ratio of rest-frame to observed frequency vq/v and how it varies along an orbit, especially around pericenter. It is 
not necessary to measure vq separately, it can be treated as a parameter to be inferred. As a result, spectroscopic 
redshifts and pulsar timings can be treated in a unified way. We argue that redshift variation over one or a few orbits, 
with special attention given to pericenter passage, provides a possible route to testing relativity using Galactic-center 
stars, or (if discovered) pulsars on similar orbits. There are two reasons why a different strategy is called for here 
than in the binary pulsar case. First, pericenter speeds of S stars are typically much higher than those of binary 
pulsars, and second, the orbital periods are too long for cumulative effects to build up. 

The observable redshift contains several different relativistic contributions affecting the orbit of the star or 
pulsar and the light traveling to the observer. For calculations, we use a four-dimensional perturbative Hamiltonian 
formulation derived from the Kerr metric. Each relativistic effect appears conveniently as a Hamiltonian term that 
can be toggled on and off to examine its import. Computation of redshift curves using this Hamiltonian approach 



was demonstrated in a previous paper (jAngelil fc Sahal 120101 ). In this paper we have developed a pipeline for 



solving the inverse problem of inferring the orbital parameters from a redshift curve. We then compute the redshift 
resolution needed to distinguish between redshift curves with and without each relativistic perturbation included, 
thus simulating the analysis of future observations. The redshift resolution needed to uncover each effect is a few 
times finer than the maximum contribution of that effect. 

In our treatment, we have assumed that frequency data is the only observable at hand. The prospects for 
detection would be improved if astrometric information were included in the fitting procedure. Information provided 
by astrometry is particularly potent in constraining the angular Keplerian elements, and would alleviate the burden 
placed on spectroscopy or pulse timing. 

A further way in which the detection prospects could be improved would be to obtain data from multiple 
orbits. The analysis of frequency data from stars or pulsars with periods < 1 yr would not only benefit from the 
steep period-dependence of the relativistic effects, but would also boost our chances of resolving cumulative effects. 
Precession effects arc naturally of this type. 

Figure|6]illustratcs the redshift accuracies required for S2. With current instrumentation, capable of ^ 10 km s^^ 
accuracies, detecting time dilation appears comfortably feasible. Detecting space curvature appears feasible with a 
modest improvement in redshift resolution. On the other hand, the discovery of a pulsar on an orbit comparable 
to S2 could push the effective redshift accuracy to < Im s~^ and, in principle, bring frame dragging and even 
quadrupole effects within reach, as shown in Figure |9l 

This paper, however, assumes a source of negligible mass in pure Kerr spacetime. The finite mass of S2 (being 
< IO^^A/bh) would perturb the redshift by a similar factor, via its effect on the motion of the supermassive black 
hole. This would be much smaller than the space-curvature effect, but more than the frame-dragging contribution. A 
potentially much more serious problem, however, is the Newtonian perturbations due to other mass in the Galactic 
centre region. This Newtonian "foreground" is not necessarily fatal - the very specific time-depe ndence of the 



relativistic effects may enable them to be extracted from under a larger Newtonian perturbation (jMerritt et al, 
20101) - but further research is needed to assess this. 



If sources inwards of S2 are discovered — someth ing we may hope for from the E-ELT and the Square Kilometre 



Array (SKA) — the prospects for relativity improve. iPfahl fc Loebl (J2004I ) argue th at there may b e 100 pulsars with 



orbital periods less than 10 years, although the number may be much smaller (e.g. iMerrittI 120091 ). A large fraction 



-10- 

of these are expected to be found by the SKA. Not only do the rclativistic effects get stronger as we move closer to 
the black hole, but the Newtonian perturbations arc likely to weaken. Figures [7] and [8] shows the orbital periods at 
which one would need to find stars or pulsars given redshift accuracies of 10 km s~^ and 1km s~^ respectively. If 
stars with a period of less than one year are discovered, the prospects become very exciting indeed. 
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Effect 


Orbit 


Light path 


Classical 


p-l/3 




Time dilation 


p-2/3 




Schwarzschild 


p-1 


P-I 


Frame-dragging 


p-3/4 


p-5/3 


Quadrupole 


p-5/3 


p-5/3 



Table 1: Relativistie effects considered in this paper, and how they scale with the orbital period. 



Effect 


Required redshift accuracy 
for S2 


Required orbital period (yr) 

for given redshift accuracy 

10 km s^i Ikms^i aOcms^i 


Time dilation 


-60kms-i 


-50 


>PS2 


>PS2 


Schwarzschild 


-3kms-i 


~13 


-30 


>i^S2 


Frame dragging 


- lOms-i 


-0.5 


-0.8 


>PS2 


Quadrupole 


- 50cms-i 


-0.04 


-0.1 


-PS2 



Table 2: Summary of the observational thresholds at which different relativistie effects are exposed, assuming a 
source of negligible mass in a pure Kerr spacetime. 
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Fig. 1. — Pericenter distance (in units of GM/c^) against orbital period for a variety of systems. The known S stars, 
having smaUer pericenter distances, are more relativistic than known binary pulsars. But the long orbital periods 
of the S stars render it infcasible to measure cumulat i ve eff ects over many orbits. Hence other tech niques must be 



devised . The pulsar examples are taken from lLorimeiJ (|2008() . For the S stars, the orbital elements in lGillessen et al. 



(|2009al ) have been used. In this paper we also treat fictitious stars that lie along the dashed line, that is, having a 
range of semi- major axis a values but with the same eccentricity and angular elements as S2. 
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Fig. 2. — Depiction of the system we deal with, shown at four shces of coordinate time. The curve ending in a 
star symbol traces a stellar orbit, starting from apocentcr, through its second pericenter passage. The star emits 
photons in all directions at equal intervals of proper time. Dots show photons which will reach the observer (marked 
x). Note that the dots represent a sequence of photons emitted at different places, hence joining the dots does not 
represent the path of any particular photon. Time dilation and Schwarzschild perturbations are included in this 
example, but not frame dragging or quadrupole. This means that the star's orbit precesses, and the photons are 
lensed. In order to make the relativistic effects visually discernible, the star has a very low value of a = 100. Such 
a star may, however, yet be discoverable by the E-ELT (JLvubenova fc Kissler-Patiai200£^ ). Also, the observer has 
been placed at an unrealistically close distance of 300. 
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Fig. 3. — Redshift contributions of different orbital effects. Shown here is the contribution of each term in Equation 
(p|) at pericenter, with the orbital period being varied and the other orbital parameters fixed at the S2 values. The 
pericenter distance varies from 3600 (approximately the value for S2) down to 6. Note that for the last two signals - 
the frame-dragging and the quadropole - we have taken the black hole spin to be maximal and to point perpendicular 
to the line of sight. The signal scalings are listed in Table [T] The curves begin to intersect as we move towards the 
strong field regime - attributable to the the breakdown in our perturbative approximation for small r. 
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Fig. 4. — Like Figure[3]but for photon propagation delays (Equation |4l). The Schwarzschild propagation signal is for 
the most part slightly smaller than the Schwarzschild orbital signal, and scales in the same way. The frame-dragging 
propagation signal is considerably smaller than the corresponding orbital signal, and scales like z^^ oc P~^/^, as 
opposed to zpD oc P^'^'^ . The frame-dragging signal remains approximately an order of magnitude weaker on the 
photons than on the star. This suggests that in attempting to measure the spin of the black hole in the post- 
Newtonian regime, its manifestation on the orbit is what matters. Note however, that this is not the case for the 
Schwarzschild effects - for which neither the photon nor orbit perturbations may be neglected. 
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Fig. 5. — Relativistic redshift effects, with orbital and light-path contributions summed. Two estimates for the 
signal due to the extended mass distribution are also shown, using the cru de model (1161) normali zed so that the 
circular velocity at the r = 10^ GMga/c^ (or -- 0.1 pc) is -- 100 km s"^ (cf. ICillessen et al.l boOQal) . The flat and 
sloping dashed curves correspond to 7 = 2.5 and 1.5 respectively. 
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Fig. 6. — The signal-to-noise ratio (as defined by equation \T7} for different relativistic effects in the orbit of S2 
as a function of redshift accuracy. Each point on this figure corresponds to a simulated data set of 200 redshifts, 
distributed over the orbit but with highest density around periccnter. 
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Fig. 7. — Similar to Figure [H except that the redshift accuracy is fixed at 10 km s~^ and the orbits vary, being 
scaled-down and speeded-up versions of the orbit S2. For large S/N these curves scale with the powers in Table [TJ 
Current instrumentation, operating under optimum conditions, should manage an accuracy of 10 km s~^ indicating 
that gravitational time dilation should be able to be detected on some of the currently known S Stars. 
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Fig. 8. — Like Figure [71 but fo r a rcdshift accuracy of 1 km s ^ 
( Lvubenova fc Kissler-Patijboogi ). 
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Fig. 9. — Like Figures [7] and [8] but for redshift accuracy 30 cm s ^. Pulse timing accuracies at this level are already 
available, known pulsars orbiting the black hole on suitably short orbits however are not. 



